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The  infrared  properties  of  thin  Pt/Al203  granular  metal-insulator  composite  films  (cermets) 
were  investigated  using  Fourier  transform  spectroscopy.  The  films  were  grown  by  dual  electron  beam 
evaporation  onto  polished  single  crystal  sapphire  substrates.  They  were  harrowed  from  Dr.  J.V. 
Mantese  of  General  Motors  Research  Labs.  The  metallic  volume  fraction,  which  was  controlled  by  the 
deposition  conditions,  ranged  from  23%  to  100%  in  the  set  of  samples.  The  percolation  threshold 
estimate  from  the  temperature  dependence  of  the  dc  resistance  is  greater  than  50%,  which  implies  that 
the  Pt  particles  are  highly  correlated.  The  percolation  threshold  for  a  random  metal-insulator  mixture  is 
less  than  20%. 
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The  transmission  and  relative  reflection  of  the  films  were  measured  at  room  temperature  using  a 
Nicolet  740  Fourier  transform  infrared  spectrometer  (FTIR),  which  arrived  at  the  beginning  of  the 
year.  The  mid  and  near  infrared  regions  (400-15250  cm-J)  were  covered  using  three  combinations  of 
sources,  beamsplitters,  and  detectors.  The  reflection  standard  was  a  100%  Pt  film  from  the  sample 
set.  The  angle  of  incidence  for  the  reflection  measurements  was  5°. 
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The  data  were  compared  with  three  effective  medium  theories  which  describe  distinct 
micros tructural  topologies:  Maxwell-Gamett  (dielectric-coated  metal  spheres),  Bruggeman  (metal  and 
insulator  treated  on  the  same  footing),  and  the  Ping  Sheng  model  (metal-coated  dielectric  and 
dielectric-coated  metal  spheres  treated  symmetrically).  The  Sheng  model  is  in  a  sense  a  combination  of 
the  other  two.  It  was  developed  specifically  for  application  to  cermet  films.  The  underlying  argument 
is  based  on  a  probabilistic  growth  mechanism  for  the  films.  The  best  available  data  on  bulk  optical 
properties  were  used  to  model  the  complex  dielectric  functions  of  the  substrate  and  the  film 
components. 


The  measured  relative  reflection  for  some  of  the  samples  is  compared  with  the  predictions  of 
the  three  models  in  the  accompanying  preprint.  Figures  1  and  2  show  a  comparison  between  the 
measured  transmission  and  the  predictions  of  the  three  models.  Samples  with  greater  than  50%  Pt  do 
not  transmit  well.  Figure  3  shows  the  volume  fraction  dependence  of  the  transmission  at  two  selected 
frequencies.  The  Sheng  model  provides  the  best  agreement  with  the  data  over  the  complete  range  of 
metallic  composition,  in  agreement  with  the  analysis  of  the  relative  reflection.  Thus,  the  spectroscopic 
measurements  suggest  that  the  Pt  particles  in  the  films  are  correlated  rather  than  randomly  dispersed. 
No  direct  corroboration  by  a  technique  such  as  electron  microscopy  has  been  obtained. 

The  most  striking  disagreement  between  theory  and  experiment  occurs  near  the  dc  percolation 
threshold,  fc.  It  is  best  observed  in  the  reflection  data.  This  result  is  not  surprising  as  effective 
medium  theories  are  known  to  fail  near  percolation.  Details  of  cluster  morphology  not  included  in  the 
effective  medium  treatments,  such  as  nonspherical  shape,  are  most  important  near  fc. 

B.  Far  Infrared  Properties  of  Small  Bismuth  Particles 


Small  bismuth  particles  were  prepared  by  inert  gas  evaporation  and  dispersed  in  paraffin, 
which  does  not  absorb  strongly  in  the  frequency  range  of  interest.  The  far  infrared  absorption 
associated  with  the  plasma  sphere  resonance  in  1  ^m  diameter  bismuth  particles  was  measured, 
identified,  and  compared  with  a  model  we  had  developed  earlier.  The  measured  absorption  coefficient, 
which  peaks  near  170  cm-l,  is  enhanced  by  about  a  factor  of  five  relative  to  the  prediction. 

The  dependence  of  the  plasma  resonance  absorption  on  particle  size  was  studied.  The 
resonance  is  not  observed  for  samples  with  mean  particle  diameters  less  than  .25  jtm.  Enhanced 
background  absorption  (i.e.,  weak  frequency  dependence)  was  observed  for  the  smallest  particles. 
Particle  sizes  were  measured  by  transmission  electron  microscopy. 

The  far  infrared  transmission  of  a  powder  of  1  nm  diameter  free  standing  bismuth  particles  in  a 
magnetic  field  was  studied.  A  group  of  cyclotron-like  resonances  (resonance  frequency  approaches 
zero  as  the  magnetic  field  is  removed)  was  observed.  The  observed  resonances  confirmed  earlier 
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Figure  2.  Transmission  spectra  of  Pt-Al203  granular  composite  films.  The  predictions  of 
the  Maxwell-Gamett  (MG),  Bruggeman  (BR)  and  Ping  Sheng  probabilistic  growth  (PG) 
models  are  also  shown.  The  PG  model  describes  the  data  best. 
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Figure  3.  Volume  fraction  dependence  of  the  infrared  transmission  of  Pt-Al2C>3  granular 
composite  films  at  selected  frequencies.  Also  shown  are  the  predictions  of  the  Maxwell- 
Gamett  (MG),  Bruggeman  (BR)  and  Ping  Sheng  probabilistic  growth  (PG)  models.  The 
PG  model  best  describes  the  data.  The  samples  do  not  transmit  above  the  percolation 
threshold,  so  the  comparison  shown  here  is  only  meaningful  for  relatively  low  metallic 
volume  fractions. 


unpublished  experimental  results  of  Chin  and  Sievers.  In  addition,  we  observed  evidence  for  many 
more  resonances  in  qualitative  agreement  with  our  model.  Better  data  (higher  resolution,  improved 
sensitivity)  are  required  to  make  the  comparison  quantitative.  We  obtained  evidence  that  the  plasma 
sphere  resonance  is  blue-shifted  in  a  magnetic  field,  in  agreement  with  theory.  We  performed  the 
magnetooptical  studies  at  Cornell  University,  where  a  cryostat  with  a  superconducting  solenoid  is 
available. 

Further  details  regarding  the  measurements  on  bismuth  particles  are  given  in  the  accompanying 
preprint  Many  other  details  remain  to  be  written  up. 

Our  expressions  for  the  extinction  coefficient  of  a  gyrotropic  sphere  in  the  long  wavelength 
limit  were  applied  to  small  bismuth  particles  in  zero  magnetic  field.  For  the  anisotropic  diejectric 
tensor  of  bismuth,  the  expressions  are  not  sufficiently  general  to  treat  the  case  of  nonzero  applied  field. 
However,  the  average  over  particle  orientations  necessary  to  treat  our  samples  is  straightforward  to 
perform. 

Improvements  in  the  far  infrared  system  over  the  past  year  include  the  construction  of  an  output 
coupler  which  uses  mirrors  to  replace  the  TPX  lens  that  came  with  the  SPECAC  interferometer.  The 
intensity  of  modulated  radiation  emitted  from  the  interferometer  was  increased,  especially  for 
frequencies  greater  than  200  cm-l,  a  region  of  interest  for  bismuth.  A  new  sample  rotator  insen, 
which  permits  multiple  exchanges  of  samples  with  minimal  accumulation  of  frozen  air  on  the  sapphire 
window  which  isolates  the  detector  chamber  from  the  rest  of  the  cryostat,  was  constructed.  A  small 
dewar  to  house  one  of  the  composite  silicon  bolometers  was  purchased  from  Infrared  Laboratories 
using  startup  money  from  the  University.  This  dewar  is  convenient  for  diagnostics  as  well  as  room 
temperature  transmission  measurements.  In  addition,  it  is  compatible  with  the  Nicolet  740,  for  which 
it  can  be  used  to  improve  performance  in  the  far  infrared. 

C.  Theory 

1.  Magnetic  Field-Induced  Resonances  in  Small  Spherical  n-InSb  Particles 

As  pan  of  a  continuing  investigation  of  the  interaction  of  electromagnetic  waves  with  a  small 
gyrotropic  sphere,  a  simple  model  for  magnetic  field-induced  resonances  in  an  isotropic  single-carrier 
polar  semiconductor  sphere  was  developed  and  applied  to  n-InSb.  Free  carriers  were  treated  using  the 
Drude  model,  which  is  characterized  by  the  plasma  frequency  cop.  The  optic  phonons  were  modeled 
by  a  Lorentz  oscillator,  for  which  the  transverse  and  longitudinal  frequencies  are  denoted  cor  and  col, 
respectively.  Depending  on  the  doping  and  temperature,  one  can  have  cop  <  cor,  cox  <  o>p  <  col,  or  cop 
>  col-  Both  electric  dipole  and  magnetic  dipole-electric  quadrupole  field  dependent  resonances  are 
predicted  in  the  long  wavelength  limit.  To  examine  magnetic  dipole -electric  quadrupole  resonances. 


use  is  made  of  our  recendy  derived  general  expressions  for  the  gyrotropic  sphere  (See  reprint).  The 
behavior  of  the  electric  dipole  resonances  is  similar  to  the  well-known  plasma  shifted  cyclotron 
resonances  with  additional  coupling  to  the  optic  phonons  (Figure  4). 

2.  The  Effect  of  Clustering  on  the  Far  Infrared  Absorption  of  Small  Metal  Particles:  A 
Simple  Model 

Clustering  of  the  particles  plays  an  important  role  in  the  anomalous  enhancement  of  the  far 
infrared  absorption  coefficient  of  small  metal  particles,  as  was  shown  experimentally  by  several 
groups1'3.  Curtin  and  Ashcroft4  introduced  a  number  of  models  which  show  how  clusters  of  particles 
produce  enhanced  absorption.  One  of  these,  the  cluster  percolation  model,  applies  to  clusters  of 
electrically  isolated  particles.  The  low  conductivity  of  a  cluster  near  the  dc  percolation  threshold  leads 
to  the  enhancement.  The  enhancement  can  be  regarded  as  a  resonance  in  the  metallic  volume  fraction 
of  the  cluster,  in  analogy  with  the  plasma  sphere  resonance  in  the  frequency  dependence  of  the  optical 
properties  an  isolated  small  particle.  Curtin  and  Ashcroft  modeled  their  clusters  using  an  approach 
based  on  the  real  space  renormaliztion  group.  Here  we  examine  a  simpler  model  which  produces  a 
similar  effect 

A  cluster  of  particles  is  replaced  by  an  homogeneous  sphere  made  up  of  an  effective  medium. 
We  choose  to  use  the  Bruggeman  model  for  the  effective  dielectric  function.  This  model  has  a 
percolation  threshold,  which  is  required  to  produce  the  enhancement.  An  actual  sample  will  contain 
clusters  with  a  range  of  metallic  volume  fractions,  with  perhaps  only  a  small  fraction  near  percolation. 
Like  Curtin  and  Ashcroft,  we  assume  that  the  details  of  the  cluster  distribution  are  unimportant  and 
consider  a  distribution  which  is  constant  over  a  range  of  volume  fractions  that  includes  the  percolation 
threshold.  Tne  volume  fraction  of  metal  in  the  sample  is  presumed  to  be  much  less  than  unity.  We  use 
the  Maxwell-Gamett  model  in  the  dilute  limit  to  obtain  the  effective  dielectric  function  for  the  medium. 
In  summary,  effective  medium  theory  has  been  applied  twice:  first  to  treat  the  clusters,  then  to  treat  the 
entire  sample. 

For  dispersed  particles  the  absorption  coefficient  is  given  by 


a  Q  =  18x  fe0  to 


e2(co) 


e,(co)  +  2eJ  +e2!toT 
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Figure  fc.  Magnetic  field  dependence  of  electric  dipole  resonances  in  small  spherical  n-InSb 
particles.  The  optic  modes  of  InSb  are  included  in  the  model.  The  solid  lines  indicate  the 
"cyclotron  resonance  active"  resonances,  the  dashed  lines  the  "cyclotron  inactive" 
resonances.  The  parameters  used  in  the  calculations  were  col  =  196  cm-1,  coy  =  184  cm1, 
6l  =  16,  and  m*/m  =  0.019.  The  plasma  frequency,  which  is  determined  by  the 

carrier  concentration  is:  a)  top  =  50  cm1,  b)  (Op  =  190  cm-',  and  c)  o)p  =300  cm-'.  The 
interaction  between  the  plasma  shifted  cyclotron  resonance  and  the  optic  phonons  is 
significant  near  2(X)  cm  ' 
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Here  f  is  the  volume  fraction  of  metal,  ei  and  e2  are  the  real  and  imaginary  parts  of  the  complex 

/v 

dielectric  function  of  the  metal,  £0  is  the  dielectric  constant  of  the  host,  and  co  is  the  frequency  in 
wavenumbers.  For  clusters. 


a(to)  =  18j:  f  e^g) 


J 


dp  P(p)  e2(gj;p) 


£1(5;p)  +  2e(]  +£2(co;p)" 


where  P(p)  is  the  distribution  of  metallic  volume  fractions  within  the  clusters  andejc^  is  the  real 
(imaginary)  part  of  the  complex  effective  dielectric  function  for  a  cluster. 

The  model  is  applied  to  100  A  diameter  silver  particles  in  KC1.  The  parameters  of  the  Drude 
model  for  silver  are  (Op  =  72,600  cm-1  (plasma  frequency)  and  1/t  =  370  cm-1  (electron  scattering  rate 
due  to  the  boundaries  of  the  particles).  The  low  frequency  dielectric  constant  of  KC1  is  eo  =  4.2. 
Figure  5  shows  the  predicted  frequency  dependence  of  the  far  infrared  absorption  coefficient  for 
dispersed  and  clustered  particles.  In  this  case,  clustering  produces  an  enhancement  of  about  a  factor  of 
ten. 

III.  RESEARCH  PLANS  FOR  1989 

A.  Infrared  Properties  of  Pt/Al203  Cermet  Films  and  Pt/KBr  Pressed  Pellets 

We  have  started  a  study  of  the  infrared  reflection  of  pressed  pellets  of  Pt  particles  (purchased  Pt 
black)  in  KBr.  A  similar  study5  was  performed  on  inert  gas  evaporated  silver  particles  imbedded  in 
KC1  by  a  group  at  Ohio  State  several  years  ago.  We  choose  to  study  Pt  particles  so  that  our  results 
can  be  compared  with  our  work  on  the  Pt-Al203  cermet  films.  The  two  systems  probably  have  very 
different  microstructures.  The  particles  in  a  cermet  film  are  apparently  correlated.  In  a  pressed  pellet, 
the  particles  should  be  randomly  dispersed  on  some  length  scale.  There  may  be  clustering  if  repeated 
grinding  and  pressing  does  not  produce  an  homogenous  mixture.  Cummings  et  al.5  concluded  that  the 
Bruggeman  model  provides  a  good  description  of  their  data.  It  will  be  interesting  to  learn  whether  the 
percolation  threshold  of  the  pressed  pellet  system  is  below  20%  metallic  volume  fraction,  in  agreement 
with  theory  for  random  mixtures  and  in  contrast  with  the  result  for  the  cermet  system. 

In  addition,  there  is  a  possibility  of  obtaining  a  new  type  of  metal-insulator  thin  film  composite 
system  in  the  coming  year  through  my  collaboration  with  Dr.  J.V.  Mantese  of  General  Motors 
Research  Labs. 
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Figure5  .  Effect  of  clustering  on  the  far  infrared  absorption  coefficient  of  100  A  Ag  particles 
embedded  in  KC1.  The  calculation  is  described  in  the  text.  For  this  simple  model,  an 
enhancement  of  about  a  factor  of  ten  is  predicted  if  the  particles  are  clumped,  a)  linear 
plot  b)  log-log  plot. 


Theoretically,  we  have  some  interest  in  investigating  the  Sheng  model.  During  a  conversation 
at  the  Fall  MRS  Meeting,  Ping  Sheng  pointed  out  that  his  model  predicts  an  infrared  plasma  edge  in  the 
reflectivity  of  cermets  that  does  not  come  out  of  the  other  oft-used  effective  medium  theories.  Will  this 
effect  be  observable  in  composite  systems  based  on  Pt  metal,  which  is  non-Drude?  For  a  Drude  metal 
and  a  transparent  host,  the  Sheng  model  predicts  a  great  detail  of  structure  in,  e.  g.,  the  reflectivity  that 
does  not  come  out  of  the  Maxwell-Gamett  or  Bruggeman  models. 

The  possibility  of  applying  new  theoretical  developments  in  the  optical  properties  of  metal- 
insulator  composite  materials,  especially  in  the  region  near  percolation,  to  the  Pt-Al203  system  will  be 
investigated.  Several  interesing  theoretical  results  based  on  scaling  theories,  the  renormalization 
group,  etc.,  were  presented  at  the  ETOPIM2  Conference  in  Paris  last  summer. 

We  plan  to  improve  our  measurements  of  the  reflectivity  of  composite  materials  by  measuring 
the  absolute  rather  than  relative  reflectivity.  We  purchased  a  V-W  reflectometer  last  year,  but  it  has  not 
been  able  to  do  the  job  thus  far.  We  have  continuing  discussions  on  this  device  with  Harrick  Scientific 
Corp.  and  hope  that  the  problems  are  resolved  in  the  next  few  months.  Once  we  can  measure  the 
absolute  reflectivity,  we  plan  to  develop  a  computer  program  to  convert  transmission  and  reflection 
measurements  on  the  thin  cermet  films  to  the  optical  constants  n  and  k. 

B .  Far  Infrared  Properties  of  Small  Bismuth  Particles 

Our  proposal  to  the  Defense-University  Research  Instrumentation  Program  (DURIP)  has 
recently  been  selected  for  funding.  A  cyrostat  with  a  superconducting  magnet  for  far  infrared 
magnetooptical  studies  and  a  far  infrared  laser  will  be  purchased  and  installed  this  year.  The  magnetic 
field  dependence  of  far  infrared  resonances  in  bismuth  particles  will  be  investigated  and  compared  with 
the  predictions  of  our  model.  An  in-house  facility  for  magnetooptics  provides  the  opportunity  for  a 
thorough  magnetooptical  investigation  of  bismuth  particles  that  was  not  possible  when  we  were 
required  to  make  short  trips  to  Cornell. 

Studies  of  the  size  dependence  of  the  far  infrared  absorption  in  bismuth  particles  will  continue. 
Besides  the  plasma  sphere  resonance,  a  new  probe,  an  interband  transition  at  higher  frequencies,  will 
become  available  if  we  are  successful  in  developing  a  low  temperature  silicon  window  and  if  the 
bandwidth  of  the  SPECAC  interferometer  is  found  to  extend  to  sufficiently  high  frequencies.  Csl  will 
be  examined  as  a  potentially  superior  host  material  over  paraffin  at  high  frequencies.  In  coversations  at 
the  ETOPIM2  Conference  in  Paris,  I  learned  that  doping  of  bulk  bismuth  has  been  used  to  modify  the 
concentration  of  carriers.  If  doped  small  particles  of  bismuth  can  be  produced,  the  doping  should 
strongly  affect  the  far  infrared  resonances,  both  in  a  magnetic  field  and  in  zero  field. 

Improvements  to  the  far  infrared  system  will  continue  as  the  need  arises. 
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Extinction  of  electromagnetic  waves  by  a  small  gyrotropic  sphere 
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The  general  solution  of  Ford  and  Werner  for  the  scattering  and  absorption  of  a  plane  electromag¬ 
netic  wave  by  a  gyrotropic  sphere  is  analyzed  in  the  long-wavelength  limit.  An  expression  is  ob¬ 
tained  for  the  extinction  cross  section  which  includes  the  electric  dipole  term  as  well  as  the  magnet¬ 
ic  dipole  and  electric  quadrupole  terms,  which  are  coupled.  The  results  agree  with  the  Mie  solution 
in  the  limit  of  vanishing  gyrotropy.  The  expressions  are  applied  to  a  lossless  single-component  plas¬ 
ma  sphere.  Agreement  is  obtained  with  recent  work  based  on  a  quasistatic  approximation. 


I.  INTRODUCTION 

Recently  Ford  and  Werner1  presented  a  solution  to  the 
problem  of  the  scattering  and  absorption  of  a  plane  elec¬ 
tromagnetic  wave  by  a  gyrotropic  sphere  made  up  of  ma¬ 
terial  having  a  frequency-dependent  complex  dielectric 
tensor  of  the  form 


~e*y  0  •  (I) 

0  0  ea 

Gyrotropy  is  frequently  introduced  via  a  static  uniform 
applied  magnetic  field  H.  Although  not  completely  gen¬ 
eral,  this  tensor  applies  to  microwave  and  far-infrared 
measurements  on  narrow-gap  semiconductors2-9  as  well 
as  electron-hole  droplets  in  Ge  (Refs.  10-23)  for  certain 
orientations  of  the  crystal  with  respect  to  the  applied 
magnetic  field. 

Approximate  solutions  which  apply  to  particles  that 
are  much  smaller  than  the  wavelength  of  the  electromag¬ 
netic  radiation  are  important  because  they  frequently  ap¬ 
ply  to  experiments  and  because  the  expressions  obtained 
in  this  limit  are  relatively  simple  and  accessible  to  physi¬ 
cal  interpretation.  Consider  the  Mie-Debye  problem24,25 
of  scattering  and  absorption  by  an  isolated  homogeneous 
isotropic  sphere  (H=0)  having  radius  a  and  complex 
dielectric  function  e(a>)  embedded  in  a  nonabsorbing  host 
with  dielectric  constant  e2 ■  Define  the  dimensionless  pa¬ 
rameter 

x  =ka  =e\n(oa  /c  (2) 


where  w  is  the  frequency  of  the  wave  and  k  is  the  wave 
number.  In  the  limit  x  «  1  the  pair  of  infinite  indepen¬ 
dent  multipole  series  reduces  to  three  terms,  to  order  x3. 
For  example,  the  extinction  cross  section  is26,27 


t  —  rrie) 

+<*\  +^2 


where 


(3) 


are  the  electric  dipole,  magnetic  dipole,  and  electric 
quadrupole  terms,  respectively.  For  most  situations  the 
electric  dipole  term  is  the  most  important,  and  the  mag¬ 
netic  dipole  and  electric  quadrupole  terms  are  the  leading 
corrections. 

In  the  long-wavelength  limit,  it  is  also  possible  to  ob¬ 
tain  expressions  for  the  extinction  coefficient  and  related 
quantities  based  on  a  quasistatic  approximation.  One  can 
distinguish  the  condition  x  «1,  i.e.,  the  wavelength  of 
the  incident  wave  in  the  host  is  much  larger  than  the  par¬ 
ticle  size,  from  the  usually  more  restrictive  condition  that 
the  wavelength  inside  the  particle  also  be  much  greater 
than  the  size  (the  Rayleigh  limit).  The  electromagnetic 
wave  is  approximated  by  uniform  oscillating  electric  and 
magnetic  fields.  Retardation  effects  are  ignored.  The 
responses  of  the  particle  to  the  electric  and  magnetic 
fields  are  treated  separately.  The  results  agree  with  the 
electric  and  magnetic  dipole  terms  of  the  Mie  series  to 
leading  order  in  x.  It  is  straightforward  to  generalize  the 
quasistatic  approximation  for  the  electric  dipole  term  to 
particles  having  ellipsoidal  shape28  and/or  anisotropic 
dielectric  tensor.29  The  quasistatic  approximation  also 
underlies  the  arguments  leading  to  effective-medium 
theories30, 31  for  composite  materials. 

Ford  and  Werner'  obtained  expressions  for  the  electric 
and  magnetic  dipole  moments  for  a  gyrotropic  sphere  in 
the  Rayleigh  limit.  Their  expressions  agree  with  results 
derived  using  a  quasistatic  approximation,32-34  for  which 
the  electric  and  magnetic  fields  are  assumed  to  decouple. 

Furdyna  et  al.iS  and  Goettig  and  Trzeciakowski36  re¬ 
cently  calculated  the  electromagnetic  modes  of  oscillation 
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of  a  small,  lossless,  single-carrier,  conducting  sphere  in  a 
uniform  static  magnetic  field.  Their  quasistatic  approxi¬ 
mation  allowed  for  coupling  of  the  electric  and  magnetic 
modes  through  displacements  of  the  plasma  with  respect 
to  the  uniform  compensating  background.  The  modes 
possess  nonvanishing  surface  and,  in  some  cases,  volume 
charge  densities.  The  interaction  of  external  electromag¬ 
netic  radiation  with  these  modes  and  the  absorbed  power 
were  determined  using  a  quantization  scheme  originally 
developed  for  the  zero-field  case.37,38  The  frequencies 
and  power  absorption  for  the  electric  dipole  resonances 
are  in  agreement  with  earlier  work,32-34  but  not  for  the 
magnetic  dipole  modes.  The  origin  of  the  discrepancy  is 
that  the  electric  and  magnetic  fields  do  not  decouple,  con¬ 
trary  to  the  assumptions  of  the  earlier  work.  In  the 
long-wavelength  limit,  the  magnetic  dipole  and  electric 
quadrupole  terms,  which  are  of  the  same  order  in  x,  mix. 
This  mixing  was  also  recognized  in  some  of  the  earlier 
work.1'7 

The  question  remains  as  to  why  the  long-wavelength 
limit  of  the  Ford-Werner  (FW)  solution  agrees  with  the 
older  incorrect  quasistatic  approximation  rather  than  the 
recent  results  of  Furdyna  et  al.iS  and  Goettig  and 
Trzeciakowski.36  In  this  paper,  we  obtain  the  correct  ex¬ 
pressions  from  the  FW  theory  for  the  magnetic  dipole- 
electric  quadrupole  extinction  cross  sections  in  the  long- 
wavelength  limit.  In  effect  we  are  generalizing  Eqs.  (4b) 
and  (4c)  to  the  gyrotropic  case.  We  do  not  treat  correc¬ 
tions  to  the  electric  dipole  absorption  that  are  of  the  same 
order  as  the  magnetic-dipole -electric-quadrupole  terms. 
In  addition  to  resolving  the  discrepancy,  this  work  gen¬ 
eralizes  the  results  of  Furdyna  et  al. 35  and  Goettig  and 
Trzeciakowski36  to  any  dielectric  tensor  given  by  Eq.  (I). 
Multicomponent  plasmas,  phenomenological  damping, 
interband  absorption,  etc.  may  be  included  in  the  treat¬ 
ment.  The  restrictions  are  those  of  the  long-wavelength 
limit.  In  cases  of  doubt  one  should  make  use  of  the  full 
FW  theory. 

The  paper  is  organized  as  follows.  Section  II  presents 
the  derivation  of  the  extinction  coefficient  from  the  FW 
theory  in  the  long-wavelength  limit.  Although  it  is  well 
understood,  the  electric  dipole  term  is  included  for  com¬ 
pleteness.  Section  III  discusses  the  results  in  the  context 
of  previous  work  on  the  single-component  Drude  plasma. 
Section  IV  summarizes  the  conclusions. 

II.  EXTINCTION  COEFFICIENT 
IN  THE  LONG-WAVELENGTH  LIMIT 

Ford  and  Werner1  write  the  extinction  cross  section  as 
[Eq.  (2.59)  in  their  paper] 

crexl=  -Im[F(k,k )-E* /  ]  E,  | 2]  (5) 

where  F(kuk2)  is  the  vector  scattering  amplitude  [Eq. 
(FW3.49)],  k  is  a  unit  vector  pointing  in  the  direction  of 
propagation  of  the  electromagnetic  wave,  and  E,  is  the 
electric  field  vector  of  the  incident  wave.  The  asterisk 
denotes  complex  conjugation.  Equations  from  the  paper 
by  Ford  and  Werner1  are  indicated  by  the  prefix  FW. 
The  quantity  of  greatest  interest  here  is  Z”",  the  ratio  of 


A  X  A  determinants  that  appears  in  F(k,k).  For  example 
[from  Eq.  (FW3.47)], 
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where  the  expressions  for  X%°  and  Y[£a  are  given  by  Eqs. 
(FW3.35)  and  (FW3.39).  The  indices  m  and  a  are 
suppressed  in  Eq.  (6).  For  details  on  the  FW  theory, 
definitions,  and  notation  the  reader  is  referred  to  the  orig¬ 
inal  paper.1 

Ford  and  Werner1  calculated  the  electric  and  magnetic 
dipole  moments  of  the  gyrotropic  sphere  in  the  long- 
wavelength  limit.  Their  results  can  be  used  to  compute 
the  extinction  cross  section  and  related  quantities.  Their 
expression  for  the  electric  dipole  term  is  correct,  but  fur¬ 
ther  examination  of  the  magnetic  dipole  term  is  required. 
(1)  Due  to  mixing  with  the  electric  quadrupole  term,  ad¬ 
ditional  terms  of  the  same  order  must  be  considered.  (2) 
The  assumption  x  «y,  where  y=qa  with  q  a  parameter 
analogous  to  the  wave  vector  inside  the  sphere  [Eq. 
(FW3.12)],  is  not  required  and  in  fact  is  the  source  of  the 
disagreement  of  computed  resonance  frequencies  with 
Furdyna  et  al* 5  and  Goettig  and  Trzeciakowski.36  We 
obtain  the  extinction  cross  section  from  Eq.  (FW3.49). 

Consider  the  eight  terms  in  Eq.  (FW3.49)  with  1,1'  <2. 
Represent  the  indices  /,  m,  and  a  by  ( l,l',m,a ).  Then 
(1,1, m,  +  )  denotes  the  electric  dipole  term.  It  is  the 
lowest-order  term  in  x  and  thus  is  expected  to  produce 
the  largest  extinction.  The  next  term,  (l,l,m,  — ),  is  the 
magnetic  dipole  term.  Although  at  first  glance  it  appears 
to  be  of  the  same  order  as  the  electric  dipole  term,  it  is  in 
fact  of  order  x3  because  is  of  order  x2.  The  term 

(2,2, m,  —  ),  the  electric  quadrupole  term,  is  also  of  order 
x3.  The  terms  (l,2,m,  —  )  and  (2,l,m,  —  )  are  of  order  x3 
and  represent  the  interaction  between  the  electric  quad¬ 
rupole  and  magnetic  dipole  terms.  Therefore,  all  four 
terms  should  be  considered  together.  In  the  limit  of  van¬ 
ishing  gyrotropy  (e.g.,  isotropic  material  with  no  dc  ap¬ 
plied  magnetic  field),  the  cross  terms  vanish.  The  terms 
(l,2,m,  +  ),  (2,l,m,  +  ),  and  {2,2, m,  +  )  are  of  higher  or¬ 
der  in  x  and  will  be  neglected.  They  represent  the  mag¬ 
netic  quadrupole  term  and  its  coupling  to  the  higher- 
order  corrections  to  the  electric  dipole  term  in  the  long- 
wavelength  limit. 

A.  Electric  dipole  term 

The  extinction  cross  section  for  the  electric  dipole 
term,  or  equivalent  expressions  such  as  the  electric  dipole 
moment,  power  absorption,  and  absorption  coefficient, 
have  been  given  correctly  to  lowest  order  by  a  number  of 
authors.1,32  34'39'40  We  give  the  expressions  here  for 
completeness,  but  do  not  present  the  derivation,  which  is 
similar  to  that  presented  below  for  the  magnetic- 
dipole-clectric-quadrupole  terms.  The  approach  fol¬ 
lowed  by  FW  is  also  correct. 
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We  find  that 


Zm  +  _ 
11  — 


ezi  +2e2 

3 

II 

o 

e±+2e2 

m  =±1 

where  e±  =  exx  ±  i  exy .  Since 

tYr1(k)*Xk-EI][Yr,i(k)Xk-Er/!E1|2]=^-|?:-E  12 


where  Et  =  E,  /  |  E, ) , 


a\e)  =  4ira2x  Im 


+2e 


6+  +2e- 


m  =0 


m  =  ±1  . 


B.  Magnetic-dipole -electric-quadrupole  terms 

We  present  a  detailed  derivation  for  the  long- 
wavelength  limit  of  the  magnetic  dipole  term.  The 
derivations  for  the  other  terms  of  the  same  order  are 
similar.  The  results  are  summarized  in  Table  I. 

From  Eqs.  (5)  and  (FW3.49)  the  contribution  to  the  ex¬ 
tinction  cross  section  due  to  the  magnetic  dipole  term 
(1,1, m,  —  )  is 


_lm)=  16JT2  .2 


«2xlm  2  Zrr[Y7i<k)*Xk-B,] 


X[Y7,(k)xk-Bf  /  (  Bj  | 2]  , 


where  the  Yftlk)  are  vector  spherical  harmonics.34,41  ~43 
Since 

1/2 

Yn(k)=  —  i  ~  kxem  (11) 


and  k-B|  =0, 

[Y7(k)*  Xk-B,][Y7i(k)Xk-B*/ 1  B,  | 2] 

=  -^|e;-B,|2.  (12) 

Now  consider  Z™-,  which  is  defined  by  Eq.  (6).  Since 
both  x«l  and  y  « 1,  we  shall  expand  in  these  two 
small  parameters.  Unlike  FW  we  shall  not  impose  x  «y. 
W  require  the  expansions  of  the  spherical  functions:44 
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X  2 
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(13c)  A’”  ~  (X)~  — 
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where 


as  well  as  the  expansions  for  a,(x),  a2(x),  a'/Hx),  and  w  e 
a2n(x),  which  are  defined  by  Eq.  (FW2.38).  The  pro¬ 
cedure  is  to  expand  the  matrices  for  Z”~  in  x  and  y  and 
manipulate  the  first  and  second  rows  of  both  the  numera¬ 
tor  and  denominator  matrices.  The  two  matrices  can  be  and 
made  identical  to  within  factors  of  kk  (i.e.,  y)  and  thus 
cancel,  leaving  a  relatively  simple  result.  The  expansions 
of  the  row  elements  are 

l-y--^  drjk),  (14a) 

Y?-a)~-±(xl-y2)drmU) ,  (14b) 
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Further  manipulation  is  required  for  Eq.  (14c).  From  the 
eigenvalue  equation  (FW3.15)  we  obtain 


Next,  we  rearrange  the  numerator  determinant  by  modi¬ 
fying  the  first  row  according  to 


(21) 


to  obtain 


d7mlk)=- 


5 

4  — m 2 


YT~^(l-Am/Cm)d{m  .  (22) 

Similarly,  modify  the  first  row  of  the  denominator  deter¬ 
minant  according  to 


X  drJx) .  (15) 

mf+ilV 


From  the  definition  of  &fmik)  [Eq.  (FW3.21)], 
^2m^)  =  y[5(4 -m2)]'/2(?m+iiV)dfm(k) 
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to  obtain 


+  3  f 


m—4 


+  ±jmir-j \w  d2m(k)  . 

(16) 


From  the  definition  of  k  [Eq.  (FW3.12)], 
r  12 


[  1  —ikW)  = 


?  aia 


A  m  y  2 

x 


to  lowest  order. 

Factors  independent  of  k  can  be  factored  from  the  first 
row  of  each  determinant,  leaving  the  numerator  and 
denominator  determinants  equal,  so  that  they  cancel. 
Thus,  in  the  long-wavelength  limit, 


■ym-  f2_  COO  Cm 

Z"  "15  c  A  ' 


O',  Y,  and  ?  are  defined  by  Eq.  (FW3.2).  After  substitu¬ 
tion  and  algebraic  manipulation, 


Substitution  for  Am  and  Cm  leads  to  the  result 
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38 
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g2(  2g„  4-  ga ) + 2g„  ga 


m  =0 


g2(g±+ga  — 6g2)+4g±gr 


(26) 


m  =±1  . 


3c2+f±+ca 

We  now  summarize  the  results.  The  extinction  cross  section  w’  vn  have  calculated  in  the  long-wavelength  limit  is 


(4irar 
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where  the  summations  are  over  the  values  of  /,  /',  m,  and  a  discussed  previously,  and 

Y"/.(k)*Xk-B,,  a=(-lf  |Y£(k)Xk-E?/ 1  E,  | 2,  a  =  (-l)'  +  1 

— Y*,.(k)*Xk-E,,  o=(-l)/,+I  X  |Y{7(k)xk-Bf/iB1|2,  o-  =(  —  1 )' 


(27) 


(28) 


Also,  /!!=  1(  1  —  2 )( 1  —4)  •  •  • . 

Table  I  lists  the  expressions  for  Z”ff.  The  9>™<T(k)  can 
be  evaluated  using 

1/2 

3  1  (Sm'El). 
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where 

k±  =  ^2~W2(kxTiky)  , 
k0  =  kt  . 


(29b) 


(30) 


III.  DISCUSSION 

The  expressions  which  we  have  derived  correctly 
reduce  to  known  results  under  the  appropriate  condi¬ 
tions.  The  extinction  cross  section  reduces  to  the  predic¬ 
tion  of  Mie  theory  [Eqs.  (3)  and  (4)]  when  the  source  of 
gyrotropy  is  removed  and  the  complex  dielectric  function 
of  the  sphere  is  assumed  isotropic.  In  particular,  we  set 
exx  =ezz  =e,  €xy=  0.  The  cross  terms  (l,2,m,  — )  and 
(2,l,m,  — )  vanish. 

The  frequencies  of  resonances  induced  by  gyrotropy 
are  obtained  from  the  zeroes  of  the  real  parts  of  the 
denominators  of  the  expressions  in  Table  I.  For  a  given 
value  of  m,  the  denominators  for  the  magnetic  dipole, 
electric  quadrupole,  and  cross  terms  are  the  same. 
Therefore,  only  one  of  these  terms  need  be  considered  for 
the  evaluation  of  the  resonance  frequencies,  but  all  are  re¬ 
quired  to  determine  the  strength  of  the  resonances. 

To  proceed  further  it  is  convenient  to  specify  a  model 


for  the  complex  dielectric  function  of  the  sphere.  The 
simplest  model  consistent  with  Eq.  (1)  is  a  single¬ 
component  Drude  plasma  in  a  uniform  static  applied 
magnetic  field  H =Hz.  The  diagonalized  dielectric  ten¬ 
sor  is  given  by 


e±(<u)=e, 


(0(0)±Q)C) 


ga(o>)=g, 


0)0) 


(31a) 

(31b) 


where 


0)„ 


Avne1 

g,m* 


1/2 


(32) 


is  the  plasma  frequency,  n  is  the  carrier  density,  e  is  the 
charge,  m*  is  the  effective  mass,  and  e,  the  dielectric 
constant  due  to  core  polarizability.  Also,  coc  =eH /(m*c) 
is  the  cyclotron  frequency  and  55=a>-)-//r  where  r  is  the 
electronic  scattering  time.  For  a  lossless  Drude  plasma, 
r— ►oo  and  w—*o).  This  case  was  studied  by  Furdyna 
et  al. 35  and  Goettig  and  Trzeciakowski.36 

The  electric  dipole  resonances  are  given  by 


Re[ea(<u)4  2e2]=0 


(33) 


for  m=0,  which  leads  to 


<^101  =<»p 


2e2  +  ei 


(34) 


for  the  resonance  frequency.  This  frequency,  which  is  in¬ 
dependent  of  the  applied  magnetic  field,  is  the  resonance 
frequency  for  electric  dipole  plasma  motion  in  zero  field. 
The  resonances  are  labeled  according  to  the  notation  of 
Furdyna  et  al.1 5  A  resonance  is  labeled  where  r=  1 
if  the  frequency  tends  to  o),,  the  /th  zero-field  mode  of  os¬ 
cillation,  in  the  low-field  limit  and  r=2  if  the  frequency 
tends  to  zero.  For  m  =  ±  1 ,  we  set 


Re[e±((n)-f-2e2]  =  0 
to  obtain 


(35) 
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l,  ±  1, 1 : 


2 

-L 

+  tu2 

2 

2 

1/2 


(36) 


For  m  —  ±2,  resonances  are  located  using 
Re(3e2+2e±)=0  . 


(46) 


This  result  is  the  well-known  plasma-shifted  cyclotron 
resonance,32,40,45,46  which  describes  the  interaction  be¬ 
tween  the  cyclotron  resonance  and  plasmon  modes. 

The  magnetic-dipole -electric-quadrupole  terms  are 
considered  together.  For  m=  0,  resonances  are  located 
from 

Re[e2(2eiX  +€**)+ 2exxea]=0  .  (37) 


For  the  lossless  Drude  plasma,  this  condition  becomes 
(02±<occo—co 2  =  0  .  (47) 

The  frequencies  are 


-r 


2 

<uf 

■y 

2 

+  &>2 

1/2 


(48) 


The  resonance  frequencies  are  the  two  real  positive  solu¬ 
tions  of 

to4  —  (co2c+co2  +co])co2  +  [(  l+e2/2e,  )col+co2p]coj=0  (38) 

where 

co22  =  2e16)2/(3e2  +  2e,)  .  (39) 

One  of  these  modes,  which  is  written  to  order  (coc/cop)2 
as 

<u2oi  -  <02  +  {o}]  ,  (40) 

is  the  cu2oi  m°de  of  Furdyna  et  al. 35  The  other  mode, 

co~cop+\co2c/cop  ,  (41) 

is  the  bulk  magnetoplasmon  of  Goettig  and 
Trzeciakowski36  (GT).  According  to  GT  only  the  com¬ 
ponent  of  the  oscillating  magnetic  field  parallel  to  the  dc 
applied  magnetic  field  interacts  with  this  mode,  which 
arises  from  the  m  =0  magnetic  dipole  term.  However,  we 
find  that  this  mode  appears  in  the  electric  quadrupole  and 
cross  terms  as  well.  To  lowest  order  in  coc/cop,  the  power 
absorbed  due  to  the  magnetic  dipole  term  agrees  with  Eq. 
(57)  of  GT  and  the  contribution  of  the  remaining  three 
terms  is  of  higher  order,  which  resolves  the  discrepancy. 
For  m  =  ±  1 ,  the  condition  for  resonances  is 


Re(3e2-t-e±-)-e„)=0  . 


(42) 


For  the  lossless  Drude  plasma,  the  resonance  frequencies 
are  given  by  the  real  positive  roots  of 


(0}±0)cC02  —  <o\ 


.°c 

T 


=  0 


(43) 


Note  that  the  cyclotron  frequency  can  be  positive  or  neg¬ 
ative  depending  on  the  sign  of  the  carriers.  This  equation 
agrees  with  Eq.  (23)  of  Furdyna  et  al . 35  for  1=2.  The  ex¬ 
pressions  for  the  roots  are  quite  complicated.  For  the 
low-field  limit, 


2 


w2,tU~w2 


coc  5  co 
T  +  32  HT2 


(44) 


and 


**2,  sgnq,  2 


1-1 

(Uc 

2 

— - 

4 

co2 

(45) 


in  agreement  with  Furdyna  et  al. 35  sgn^  denotes  the  sign 
of  the  carriers. 


The  power  absorbed  can  also  be  calculated  for  the 
“lossless”  plasma.  It  is  related  to  the  extinction  cross 
section  by 

P=crtxiS  (49) 

where  S  is  the  magnitude  of  the  Poynting  vector  of  the 
incident  wave.  Although  our  results  do  not  appear  simi¬ 
lar  at  first  glance  they  are  identical  to  those  obtained  by 
Furdyna  et  al.,25  who  followed  a  quantum-mechanical 
approach:  (1)  Compute  the  modes  of  the  system.  (2) 
Write  down  the  Hamiltonian.  (3)  Quantize  it.  (4)  Use  it 
to  compute  the  power  absorption  by  treating  the  interac¬ 
tion  of  the  modes  with  a  perturbing  electromagnetic  field. 
Our  approach  is  classical.  We  take  our  results  for  the  ex¬ 
tinction  cross  section,  plug  in  the  dielectric  tensor  for  the 
single-component  Drude  plasma,  and  take  the  lossless 
limit  by  letting  the  relaxation  time  become  infinite.  Use 
is  made  of  the  relation 

lim  — ~ — -=ir6(jc)  .  (50) 

a-o x2+a2 

Our  results  not  only  confirm  the  work  of  Furdyna 
et  al.,25  but  are  considerably  more  general  because  we 
can  treat  a  sphere  having  any  dielectric  tensor  consistent 
with  Eq.  (1).  We  may  include  phenomenological  damp¬ 
ing,  multicomponent  systems  (which  allow  for  many  reso¬ 
nances,  including  types  not  seen  for  the  simple  example), 
optical  phonons,  and  anisotropic  carrier  effective-mass 
tensors  with  the  proper  symmetry.  The  generality  offers 
the  possibility  of  comparison  with  data  on  real  systems 
including  powdered  semiconductors  (e.g.,  n-type  InSb,  n- 
type  InAs),2-9  electron-hole  droplets  (for  specific  orienta¬ 
tions  only),10-22  and  possibly  the  characterization  of  the 
new  semiconductor  “microdots.”47-51  Cyclotron  reso¬ 
nance  and  magnetospectroscopy  in  general  has  been  of 
use  in  the  study  of  bulk  and  thin-film  semiconductors  as 
well  as  quantum  wells,  heterojunctions,  and  the  two- 
dimensional  electron  gas.52  If  quantum  dot  structures 
take  on  some  importance  for  new  device  structures,  our 
expressions  should  be  of  use  in  the  interpretation  of  mag¬ 
netooptical  data,  particularly  if  the  dots  are  spheres. 

The  behavior  of  the  electric  dipole  extinction  in  an  ap¬ 
plied  magnetic  field  in  the  long-wavelength  limit  is  well 
known.  What  is  new  is  the  expression  for  the  magnetic- 
dipole-electric-quadrupole  extinction  cross  section.  Ex¬ 
perimentally,  resonances  have  been  observed  for  electric 
dipole  excitations  in  several  systems,2-22  but  there  have 
been  few  reports  of  magnetic  dipole  excitations.20  The 
magnetic-dipole-electric-quadrupole  term  is  relatively 
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Frequency  (cm-1) 

FIG.  1.  Normalized  extinction  cross  section  for  a  small  gyro- 
tropic  sphere.  Parameters  for  a  single-component  Drude  plas¬ 
ma  are  chosen  to  model  a  1 -pm- radius  n-type  InSb  sphere  in 
vacuum:  n=5xl0,f  ein~J,  m* /m  =0.014,  l/r=0.5  cm'1, 
e,  =  1 8,  e2  =  1 .  A  6-T  magnetic  field  is  applied  in  the  z  direction. 
The  incident  wave  propagates  in  the  x  direction,  (a)  Electric  di¬ 
pole  (ED),  solid  line,  and  magnetic-dipole-electric-quadrupole 
(MDEQ),  dashed  line,  contributions,  (b)  Total  extinction  cross 
section  in  the  long-wavelength  limit.  The  ED  term  dominates, 
especially  for  the  low-frequency  (<150  cm-1)  group  of  reso¬ 
nances. 

small  and  should  be  more  difficult  to  observe,  but  the 
effect  should  be  observable.  Figure  1  shows  the  normal¬ 
ized  extinction  cross  section  aeM/(ira2),  with  parameters 
for  the  single-component  Drude  plasma  chosen  to  model 


a  1 -/am -radius  sphere  of  n-type  InSb  in  vacuum.  Figure 
1(a)  shows  the  contributions  of  the  electric  dipole  (ED) 
and  magnetic-dipole-electric-quadrupole  (MDEQ)  terms 
separately.  Figure  1(b)  shows  the  total  extinction  cross 
section  in  the  long-wavelength  limit.  The  6-T  magnetic 
field  is  applied  in  the  z  direction.  The  incident  elec¬ 
tromagnetic  wave  propagates  in  the  x  direction.  The 
MDEQ  resonances  are  weaker  by  2-3  orders  of  magni¬ 
tude  in  the  low-frequency  multiplet  below  150  cm-1,  and 
some  of  them  overlap  the  ED  resonances.  The  ED  and 
MDEQ  resonances  are  closer  in  magnitude  for  the  group 
at  high  frequencies  (>400  cm-1).  The  strength  of  the 
MDEQ  resonances  relative  to  the  ED  resonances  in¬ 
creases  with  particle  size,  but  the  applicability  of  the 
long-wavelength  approximation  becomes  more  question¬ 
able  as  well.  The  scattering  rate  was  assigned  a  low  value 
[l/r(cm-l)]=0.5  cm-1  =  l/[2ircr(s)]  to  produce  sharp 
resonances  on  the  figure.  If  the  resonances  in  actual  ma¬ 
terials  are  broader,  MDEQ  resonances  will  be  more 
difficult  to  resolve  and  identify.  For  an  experiment,  it  is 
important  to  have  well-dispersed  spherical  particles  with 
a  controlled  sharp  size  distribution.  If  the  effective-mass 
tensors  of  the  carriers  are  anisotropic,  it  is  desirable  that 
the  particles  be  aligned. 


IV.  CONCLUSIONS 

We  have  obtained  the  long-wavelength  limit  of  the  ex¬ 
tinction  cross  section  for  a  gyrotropic  sphere  from  the 
general  solution  of  Ford  and  Werner.1  We  applied  the  re¬ 
sults  to  a  lossless  single-component  Drude  plasma  sphere 
and  obtained  agreement  with  recent  calculations  by  Fur- 
dyna  et  al.is  and  GT,56  which  were  performed  using  a 
quasistatic  approximation.  Thus,  our  results  resolve  the 
discrepancy  between  these  recent  studies  and  the  earlier, 
incorrect  long-wavelength  expressions.  In  addition,  their 
generality  permits  potential  applications  to  a  number  of 
actual  materials. 
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Abstract 

We  have  observed  low  frequency  (  <  50  cm"')  magnetic  field  dependent  resonances  and  higher 
frequency  resonances  at  zero  field  in  the  far  infrared  (FIR)  absorption  of  small  bismuth  particles.  A 
model  based  on  the  quasistatic  approximation  shows  both  good  agreement  with  previous  measurements 
of  resonance  frequencies  vs.  field  and  qualitative  agreement  with  our  data.  Additional  work  is  required 
to  quantify  the  comparison. 

Introduction 

Bismuth  is  often  used  as  a  model  system  for  solid  state  plasmas  because  many  of  its  important 
energies  (Ep,  ftwp,  etc.)  are  conveniently  close  together  in  the  FIR.  The  properties  of  Bi  make  it 
especially  interesting  for  FIR  magneto-optical  studies  because  the  interaction  of  the  cyclotron-like 
resonances  with  the  plasma  modes  can  be  seen  using  currently  achieveable  fields.  The  low  carrier 
density  also  implies  a  large  Kubo  gap  [lj,  enhancing  the  possibility  of  seeing  quantum  size  effects 
(QSE)  or  metal-insulator  transitions. 

One  of  the  earliest  studies  of  small  Bi  particles  was  carried  out  by  Chin  and  Sievers  [2]. 
Motivated  by  the  idea  that  free  Bi  particles  appeared  to  align  the  bisectrix  axis  along  the  applied 
magnetic  field  [3],  they  performed  FIR  transmission  measurements  on  free-standing  Bi  particles  and 


found  three  field  dependent  resonances.  No  successful  explanation  of  these  lines  was  given.  We  have 
developed  a  model  using  bulk  Bi  parameters  that  agrees  well  with  the  data.  The  model  also  predicts 
additional  resonances  in  and  above  the  frequency  range  covered  by  Chin.  To  search  for  these 
additional  resonances,  we  investigated  the  FIR  absorption  of  small  Bi  particles  both  at  low  frequencies 
in  a  magnetic  field  and  at  higher  frequencies  at  zero  field. 

Theory 

Chin's  data  (Figure  1)  can  be  successfully  modelled  by  treating  the  Bi  particles  as  small 
spheres  with  an  anisotropic  dielectric  tensor  6.  The  quasistatic  boundary  value  problem  [4]  is  solved  to 
find  a  (the  absorption  coefficient  )  ; 

a(u)  =  §  f^  (f)  lm|t* •  B'1-  E}  +  I  fd2^(^)3  lm{B*-  F1'*-  (1) 

where  u  is  the  angular  frequency.  the  dielectric  constant  of  the  surrounding  medium,  6  t he 
dielectric  tensor  of  the  sphere,  d  the  particle  diameter,  and  f  the  volume  fraction  of  Bi.  The  tensors 
in  the  expression  are  given  by  B  =  I  +  2C0S'1  and  S_1=  Tr(8'1)l  —  (S_1)T 

To  get  8  for  a  small  Bi  sphere,  we  use  an  anisotropic  Drude  model  with  bulk  parameters  [5] 
for  the  four  carriers  : 

*  =  K  -  ^£ri>iK)'1  (2) 

j=> 

Here  n-  is  the  number  density  and  m.*  the  effective  mass  tensor  for  the  carrier.  6l  is  a  frequency 
independent  tensor  [5],  q  the  electronic  charge,  and  uf  =  w  +  £.  [Note:  The  magnetic  dipole  term  (2n 
term  in  Eq.  1)  is  incorrect  due  to  neglect  of  the  electric  quadrupole  term  [6].  However,  the  large  value 
of  SL  for  Bi  makes  the  corrections  small.) 

Figure  1  compares  the  Chin  data  with  the  theory.  The  data  agree  remarkably  well  with  three 
of  the  theoretical  curves,  allowing  us  to  identify  two  of  the  curves  as  electric  dipole  resonances,  and  one 
as  magnetic  dipole.  The  model,  however,  also  predicts  a  large  number  of  lines  that  aren’t  seen  in  these 
data  as  well  as  resonances  that  should  exist  at  higher  frequencies,  even  with  H=0. 

Experiment 

We  used  two  approaches  to  test  the  model.  First  we  looked  for  additional  low  frequency  field 
dependent  resonances  in  free-standing  powders.  The  powders  themselves  limited  the  upper  cutoff 
frequency  to  around  50  cm"'.  To  search  for  higher  frequency  resonances  we  used  pressed  pellets  of  Bi 
in  paraffin.  This  enabled  us  to  u«e  low  volume  fractions  of  Bi  and  raise  the  cut-off  frequency. 


The  Bi  particles  were  prepared  by  inert  gas  evaporation  in  Ar  [7].  Size  distributions  were 
obtained  from  transmission  electron  micrographs.  Electron  diffraction  patterns  indicated  that  the 
particles  were  crystalline. 

Samples  of  free-standing  powders  were  made  by  spreading  a  thin,  approximately  uniform  layer 
of  powder  onto  a  Mylar  sheet  and  then  sandwiching  it  with  another  layer  of  Mylar  to  keep  it  in  place 
during  sample  insertion. 

To  prepare  pressed  pellets,  the  Bi  powder  was  first  mixed  with  paraffin  powder  and  manually 
shaken.  The  paraffin  was  then  melted  and  the  suspension  of  Bi  particles  was  stirred.  After 
resolidifying,  the  sample  was  repeatedly  ground  in  a  freezer  mill  at  77*K  and  then  pressed  into  half  inch 
diameter  pellets.  The  grinding/pressing  cycle  was  repeated  3-4  times. 

The  transmission  of  the  free-standing  powder  was  measured  in  the  Faraday  geometry  for 
various  magnetic  fields  using  a  Michelson  interferometer  over  the  range  10  -  50  cm'*.  For  the  pressed 
pellets,  we  measured  the  FIR  absorption  coefficient,  again  with  a  Michelson  but  over  the  range  30  -  300 
cm'*. 

The  results  of  our  study  on  a  free-standing  Bi  powder  appear  in  Figure  2  as  ratios  of  spectra 
with  field  to  spectra  without  field.  Arrows  indicate  the  locations  of  possible  resonances,  but  many  of 
the  identifications  must  be  considered  tentative  because  of  their  small  size  (slightly  larger  than  noise  in 
the  spectra).  Even  discounting  some  of  the  smaller,  less  trustworthy  resonances,  the  data  still  show  an 
abundance  of  field  dependent  absorption  peaks  providing  at  least  qualitative  support  for  our  model. 

Quantitative  comparison  on  the  other  hand  is  difficult.  In  addition  to  the  small  size  of  many 
resonances,  interpretation  of  the  data  is  further  complicated  by  the  changing  size  of  the  peaks  as  they 
move  upward  in  frequency.  This  makes  it  difficult  to  unambiguously  trace  the  field  dependence  of  each 
resonance  from  curve  to  curve,  especially  in  regions  where  the  magnetic  field  sampling  interval  is  more 
than  lkG.  Further  confusion  results  from  the  low  transmission  level  of  the  samples  which  both 
decreases  SNR  (raising  doubts  about  the  smaller  resonances)  and  limits  the  frequency  range  over  which 
resonances  can  be  tracked.  A  direct  quantitative  evaluation  of  the  model  beyond  the  agreement  with 
Chin’s  data  must  wait  until  additional  experiments  in  progress  are  completed. 

It  is  possible  however  to  obtain  support  for  our  model  from  the  pressed  pellet  data  shown  in 
Figure  3.  The  absorption  coefficient  at  zero  field  shows  a  resonance  within  7%  of  the  predicted  location 
for  the  zero  field  combined  dipole  peak  for  randomly  oriented  particles.  The  theory  predicts  two  peaks 
but  it  is  possible  that  the  resonances  are  too  broad  to  be  resolved. 

Conclusions 

We  have  developed  a  model  for  the  far  infrared  absorption  of  small  Bi  particles  that  agrees 


with  the  data  of  Chin  [2].  Our  own  experimental  evidence  of  a  targe  number  of  field  dependent 
resonances  below  50  cm**  and  a  strong  resonance  around  180  cm'*  at  zero  field  also  suggests  at  least  a 
qualitative  agreement  with  the  model.  A  direct  quantitative  comparison  awaits  the  collection  of  more 
data. 
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Figure  Captious 


Fig.  1.  Field  dependent  resonances  for  4000A  diameter  free  standing  bismuth  particles.  The  symbols 
show  data  while  the  lines  represent  calculations  based  on  our  model.  Solid  lines  indicate  magnetic 
dipole  curves  while  dashed  lines  are  electric  dipole. 

Fig.  2.  Relative  transmission  spectra  of  5000A  diameter  free  standing  bismuth  particles  for  five 
different  values  of  the  magnetic  field.  The  dashed  lines  are  at  the  100%  level  for  each  curve.  All  the 
curves  are  plotted  with  the  same  vertical  scale.  Arrows  indicate  possible  locations  of  resonance  peaks. 

Fig.  3.  FIR  absorption  coefficient  for  1  pm  bismuth  particles  in  paraffin.  The  absorption  coefficient 
for  paraffin  is  shown  for  comparison. 
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ABSTRACT 

The  absorption  of  far  infrared  radiation  by  small  metal  particles  and  metal- insulator 
composite  materials  is  reviewed.  The  possibility  of  observing  quantum  size  effects  is 
examined.  Recent  experimental  and  theoretical  progress  in  understanding  the  anomalous 
enhancement  of  the  far  infrared  absorption  coefficient  is  discussed. 


I.  Introduction 


This  paper  briefly  reviews  the  present  understanding  of  the  interaction  of  far  infrared 
(FIR)  radiation  with  small  metal  particles  and  metal- insulator  composite  materials.  The 
important  issues  of  quantum  size  effects  (QSE),  the  anomalous  enhancement  of  the  far 
infrared  absorption  coefficient,  and  other  recent  activity  are  discussed.  This  paper  is  an 
update  to  several  more  detailed  reviews  [1-4]  of  this  topic. 


II.  Quantum  Size  Effects 

Interest  in  the  far  infrared  absorption  by  small  metal  particles  was  initially  motivated 
at  least  in  part  by  the  possibility  of  observing  transitions  between  discrete  levels.  The 
effects  of  confinement  on  the  energy  levels  of  the  conduction  electrons  in  a  particle  can 
be  described  by  quantum  box  models,  as  first  discussed  by  Frohlich  [5]  and  later  refined 
by  others  [6-10].  The  levels  calculated  in  these  models  tend  to  have  a  high  degree  of 
degeneracy  due  to  the  assumed  symmetry  of  the  particle  (usually  a  cube  or  sphere). 
Details  in  the  optical  structure  predicted  by  these  models  are  not  presently  observable 
due  to  imperfections  in  the  shapes  of  particles,  surface  roughness,  and  the  inevitable 
distribution  of  sizes. 

Kubo  [11]  introduced  a  statistical  approach  to  QSE’s  in  ultrafine  metal  particles.  He 
estimated  the  level  spacing  at  the  Fermi  energy  Ef  of  a  particle  by  inverting  the  density 
of  electron  states  at  Ef  for  a  given  spin  to  obtain  the  “Kubo  gap”  S  —  4Ep/3nV,  where 
n  is  the  density  of  free  carriers  and  V  the  volume  of  the  particle.  Using  free  electron 
values,  6  is  on  the  order  of  10cm-1  (1  meV)  for  typical  metals  (Al,  Ag,  Au,  Cu,  ...).  An 
exceptionally  large  6  is  predicted  for  Bi  particles  [12]. 

Gor’kov  and  Eliashberg  (GE)  [13]  asserted  that  Kubo’s  use  of  the  Poisson  distribution 
to  describe  the  level  statistics  neglects  correlations  between  the  levels.  According  to  GE, 
the  randomness  of  the  level  distribution  in  an  ensemble  of  metal  particles  introduced  by 
surface  roughness  is  described  by  the  statistical  theory  of  levels  (random  matrix  theory) 
originally  developed  by  Wigner  [14],  Dyson  and  Mehta  [15],  and  others  for  applications 
in  nuclear  physics.  After  corrections  for  a  number  of  errors  in  the  original  paper  [16-18], 
the  GE  model  predicts  oscillations  in  the  frequency  dependence  of  the  FIR  absorption 
coefficient  superimposed  on  a  quadratic  background.  When  account  is  taken  of  the 
particle  size  distribution  present  in  actual  samples,  the  signatures  of  QSE’s  disappear 
[18],  At  present  no  evidence  has  been  obtained  for  QSE’s  in  small  metal  particles  by  a 
far  infrared  experiment.  Thus,  theory  and  experiment  are  presently  in  agreement,  but  a 
rigorous  test  of  the  theory  requires  the  ability  to  manufacture  particles  with  a  very  sharp 
size  distribution. 

Devaty  and  Sievers  [19]  reexamined  the  possibility  of  observing  QSE’s  by  absorption 
spectroscopy  and  found  reduced  absorption  for  frequencies  below  the  mean  Kubo  gap. 
The  effect  is  weak,  but  persists  even  in  the  presence  of  a  size  distribution. 


An  anomalously  low  microwave  conductivity  recently  observed  [20]  for  small  indium 
particles  dispersed  in  oil  has  been  interpreted  as  evidence  for  a  QSE  [21]  using  the  GE 
model. 

Shklovskii  [22]  discussed  a  modification  to  GE’s  argument  based  on  electron-electron 
interactions  and  concluded  that  the  frequency  dependence  of  the  FIR  absorption  coeffi¬ 
cient  predicted  by  GE  must  be  reduced  by  a  factor  of  hw/S. 

Random  matrix  theory  has  been  an  active  area  of  research  recently  [23].  Efetov  [24] 
derived  the  two-level  correlation  functions  for  Dyson’s  three  circular  ensembles  using  the 
mathematics  of  supersymmetry.  Several  groups  have  questioned  the  assumptions  under¬ 
lying  the  GE  model  [25-28],  particularly  whether  randomness  introduced  by  boundary 
conditions  can  be  described  by  random  matrices.  The  issue  appears  to  be  unresolved  at 
this  time. 

III.  The  Anomalous  Far  Infrared  Absorption 

Most  measurements  of  the  far  infrared  absorption  coefficient  of  small  metal  particles 
imbedded  in  an  insulating  host  have  been  performed  on  samples  having  a  low  metallic 
volume  fraction  f.  In  the  limit  of  low  f,  long  wavelength,  and  assuming  that  the  particles 
are  randomly  dispersed,  the  predictions  of  effective  medium  theories  reduce  to  a  simple 
expression  which  can  also  be  obtained  from  the  Mie- Debye  solution  for  a  single  sphere. 
For  a  Drude  metal  particle  imbedded  in  a  nonabsorbing  host,  with  the  bulk  relaxation 
time  replaced  by  vp/a  (vp  =  Fermi  velocity,  a=particle  radius),  the  FIR  absorption 
coefficient  a  reduces  to  the  sum  of  electric  (ED)  and  magnetic  (MD)  dipole  terms  with 
the  properties:  1)  quadratic  frequency  dependence,  2)  size  dependence:  apD  ~  a-1  and 
<*MD  ~  a3;  3)  For  a  typical  metal,  say  Al,  the  MD  term  dominates  for  diameters  greater 
than  50A. 

An  enhancement  of  one  to  four  orders  of  magnitude  in  the  FIR  absorption  coeffi¬ 
cient  of  small  metal  particles  and  metal-insulator  composites  with  respect  to  theoretical 
predictions  was  discovered  by  Tanner  et  al.  [29]  and  subsequently  characterized  [30-33]. 
The  samples  under  study  were  typically  particles  prepared  by  inert  gas  evaporation  [34], 
sometimes  in  the  presence  of  Oj  to  produce  an  oxide  coating,  and  dispersed  in  an  al¬ 
kali  halide  host  by  repeated  grinding  and  pressing  into  pellets.  Free-standing  particles 
(smokes)  were  also  examined.  The  principal  results  are  the  ubiquity  of  the  enhancement 
(observed  for  Al,  Ag,  Au,  Cu,  Pd,  Sn,  and  Pt),  an  approximately  quadratic  frequency 
dependence,  and  a  linear  dependence  on  f  for  small  f  (with  the  exception  of  Al[33]).  In 
addition  to  the  enhancement,  the  far  infrared  properties  of  superconducting  particles 
showed  anomalous  behavior  in  the  region  of  the  gap  [35]. 

A  number  of  mechanisms  were  proposed  in  attempts  to  explain  the  anomalous  absorp¬ 
tion.  Mechanisms  intrinsic  to  isolated  particles  include  vibrations  [36-37],  eddy  currents 
and  nonlocal  effects  [38-40],  particle  size  distributions  [41],  Coulomb  effects  [42],  and 
quantum  size  effects  [43].  Mechanisms  for  enhancements  induced  by  clustering  include 


oxide-coated  clusters  [44,45j,  eddy  current  losses  [46],  and  geometrical  effects  (shape, 
close  approach  of  spheres)  [46].  These  mechanisms  do  not  explain  the  enhancement  in 
a  satisfactory  manner,  either  because  the  experimental  result  is  larger  than  theory  or 
because  the  proposed  explanation  is  too  specific  to  cover  all  the  experiments. 

An  experimental  breakthrough  occurred  with  the  realization  of  the  importance  of 
well-characterized  samples  with  controllable  properties  to  distinguish  among  proposed 
mechanisms.  Devaty  and  Sievers  [47]  developed  a  novel  composite  material,  ~  1 00 A 
Ag  particles  supported  in  gelatin,  that  could  be  sectioned  with  an  ultramicrotome  for 
examination  by  transmission  electron  microscopy  as  well  as  pressed  into  pellets  for  FIR 
measurements.  Samples  could  be  prepared  with  the  particles  well-dispersed  or  deliber¬ 
ately  agglomerated.  Studies  of  this  material  provided  evidence  for  the  important  role  of 
clustering  to  the  enhancement.  Bounds  were  placed  on  the  magnitude  of  the  enhance¬ 
ment  for  dispersed  particles  (Absorption  by  the  gelatin  dominated  absorption  by  the 
particles).  Deliberate  clustering  produced  an  increase  in  absorption. 

Additional  experiments  provided  further  important  clues.  Curtin  et  al.  [48]  showed 
that  a  heat  treatment  (melting  the  particles)  eliminated  the  anomalous  behavior  of  super¬ 
conducting  Sn  particles  near  the  gap.  This  experiment  motivated  the  models  of  Curtin 
and  Ashcroft  [49],  which  provide  an  explanation  for  the  enhancement.  Lee  et  al.  [50] 
introduced  a  new  host,  DLX-6000,  a  teflon-like  material,  which  can  be  ultramicrotomed 
for  electron  microscopy  but,  unlike  gelatin,  is  a  weak  FIR  absorber.  They  studied  fairly 
large  Ag  particles  imbedded  in  teflon  and  were  able  to  obtain  agreement  with  theory  for 
dispersed  particles  and  demonstrated  enhanced  absorption  for  a  clumped  sample. 

The  results  of  these  experiments  directed  theorists  to  focus  on  clustering  as  a  mech¬ 
anism  for  enhanced  FIR  absorption.  Curtin  and  Ashcroft  [49]  introduced  three  different 
models.  In  the  fused  cluster  model,  a  dense  cluster  is  modeled  as  a  metallic  sphere  with  a 
scattering  time  determined  by  the  radius  of  an  individual  particle.  The  short  relaxation 
time  leads  to  enhanced  eddy  current  (MD)  absorption.  The  cluster  percolation  model 
treats  a  cluster  as  a  sphere  made  up  of  an  effective  medium.  The  sample  is  modeled 
as  a  dilute  mixture  of  clusters  which  make  up  a  second  effective  medium.  Enhanced 
absorption  occurs  for  clusters  with  f  near  the  dc  percolation  threshold  fc\  i.e.,  there  is  a 
resonance  in  f.  The  enhanced  ED  absorption  is  caused  by  the  tortuous,  poorly  conducting 
clusters  near  /c.  Although  Curtin  and  Ashcroft  modeled  their  cluster  using  a  treatment 
based  on  the  real  space  renormalization  group,  any  effective  medium  theory  with  a  perco¬ 
lation  threshold  should  do  the  job,  at  least  qualitatively.  Perhaps  the  simplest  approach 
is  to  model  the  clusters  using  the  Bruggeman  model  [51]  and  use  the  Max  well- Gar  nett 
model  [52]  to  average  over  the  clusters.  The  third  model,  the  cluster-tunnel  junction 
model,  applies  to  closely  spaced  oxide-coated  particles  in  clusters.The  absorption  mecha¬ 
nism  is  photon-induced  electron  transfer.  Hui  and  Stroud  [53]  modeled  a  cluster  using  a 
self-similar  effective  medium  theory  [54]  and  examined  a  possible  role  for  tenuous  fractal 
clusters  in  the  enhanced  absorption.  Niklasson  et  ad.  [55]  also  examined  fractal  clusters. 
Claro  and  Fuchs  [56]  modeled  clusters  of  particles  by  introducing  a  distribution  of  effec¬ 
tive  depolarization  factors.  They  obtained  an  enhancement  of  3-4  orders  of  magnitude 
in  the  FIR  ED  absorption. 


In  addition  to  clustering,  recent  theoretical  work  has  considered  other  mechanisms 
including  surface  phonons  (57],  electron-phonon  coupling  [58],  diffuse  surface  scattering 
[59],  relaxation  time  effects  [60],  and  quantum  size  effects  [61]. 

Recent  experimental  work  in  the  far  infrared  has  focussed  on  new  issues.  Previous 
studies  at  low  temperature  ( T  <  20 K)  showed  no  evidence  for  temperature  dependence 
in  a.  Noh  et  al.  [62]  have  measured  the  temperature  dependence  of  a  for  a  Ag-teflon 
composite  from  room  temperature  down  to  liquid  He  temperatures  and  observed  a  small 
(~  10%)  effect.  They  found  satisfactory  agreement  with  a  model  which  included  the 
effects  of  oxide  coats,  particle  size  distribution,  and  modification  of  the  electronic  relax¬ 
ation  time  due  to  impurities  within  the  particles.  Lee  et  al.  [50,  63]  showed  that  oxide 
coatings  lead  to  enhanced  absorption  in  the  mid  and  far  infrared.  Sherriff  and  Devaty  [12] 
have  studied  the  unusual  FIR  properties  of  Bi  particles.  Kuroda  et  al.  [64]  have  measured 
FIR  absorption  by  small  particles  of  the  high- 7);  superconductor  YBa2Cu306.74-  Sievers* 
group  [65]  applied  the  Bruggeman  model  [51]  to  FIR  data  on  sintered  pellets  of  high- Tc 
superconductors  and  related  compounds. 

The  role  of  clustering  in  the  anomalous  enhancement  of  the  FIR  absorption  coefficient 
appears  to  be  understood.  However,  there  are  some  unresolved  issues.  The  dependence 
of  a  on  particle  size  has  not  been  systematically  studied.  The  inevitable  size  distributions 
make  such  a  study  difficult,  but  theoretical  predictions  should  be  tested  experimentally. 
In  addition,  there  have  been  few  FIR  studies  of  metal- insulator  composite  materials  over 
a  large  or  complete  range  of  composition  [66].  Such  studies,  which  might  focus  on  the 
region  of  dc  percolation,  would  test  effective  medium  theories  and  extend  similar  studies 
in  the  ir-vis-uv  into  a  new  spectral  region. 
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ABSTRACT 


The  reflection  of  mid-  and  near- infrared  radiation  (400- 15000 cm"  *  by  thin  Pt/Al203 
cermet  films  was  measured  using  a  Fourier  transform  spectrometer.  The  data  were 
compared  with  predictions  of  three  models  for  the  effective  optical  constants  of 
heterogeneous  materials:  Maxwell-Garnett,  Bruggeman,  and  a  simplified  version  of  a 
probabilistic  growth  model  due  to  Sheng.  Sheng’s  model  provides  the  best  description  of 
the  data  over  the  complete  range  of  metallic  volume  fraction.  This  result  is  expected 
based  on  the  most  likely  topology  of  the  films  and  is  in  agreement  with  other  work  on 
similar  systems  at  higher  frequencies. 

INTRODUCTION 


The  relationship  between  the  microstructural  topology,  the  optical  properties  of  a 
heterogeneous  system,  and  their  description  by  an  homogeneous  effective  dielectric 
function  has  been  made  clear  in  recent  years  [1].  The  prediction  of  the  optical  response 
of  a  granular  metal/insulator  composite  material  over  the  complete  range  of  composition 
provides  an  especially  stringent  challenge  to  effective  medium  theories. 

In  this  paper,  we  compare  the  predictions  of  three  simple  effective  medium  models, 
which  describe  different  microstructural  topologies,  with  the  measured  relative  reflectivity 
of  thin  Pt/Al203  films  in  the  mid-  and  near- infrared.  These  models  have  been 
successfully  applied  to  similar  cermet  systems  in  the  near  infrared  and  visible  regions 
[1,2].  We  find  that  a  simplified  version  of  a  probability  growth  model  introduced  by 
Sheng  [2,3]  provides  the  best  agreement  with  the  data  over  the  complete  range  of  volume 
fraction.  This  result  is  not  surprising  since  this  model  was  developed  specifically  for 
application  to  cermet  films.  The  greatest  disagreement  between  theory  and  experiment 
occurs  for  metallic  volume  fractions  near  the  percolation  threshold  for  the  dc 
conductivity.  Effective  medium  theories  are  known  to  fail  in  this  region[4], 

EXPERIMENTAL 


The  Pt/Al203  granular  films  were  prepared  by  coevaporation  of  Pt  and  AI2O3  onto 
single  crystal  sapphire  wafers  in  a  dual  e-beam  evaporator  with  a  low  base  pressure  (mid 
1(T’  torr).  The  volume  fraction  of  Pt  in  the  films,  f,  which  ranges  from  0.23  to  1.00, 
was  determined  by  monitoring  the  relative  deposition  rates  of  the  two  materials  using 
quartz  crystal  oscillators  inside  the  evaporation  chamber  [5].  Film  thicknesses,  measured 
with  a  Dektak  surface  profilometer,  ranged  from  1 100  to  1800  K.  The  orientations  of  the 
surfaces  of  the  sapphire  substrates  were  determined  by  back  reflection  Laue  photographs 
to  be  about  20°  from  the  ordinary  axis.  The  percolation  threshold,  fc,  as  determined  by 
the  temperature  dependence  of  the  dc  resistivity,  is  in  the  range  of  50-59%  Pt  by  volume 
[6].  Although  no  direct  determination  of  microstructure  by  electron  microscopy  has  been 
performed  on  the  films,  the  high  value  of  fc  is  a  strong  indicator  of  a  coated-grain 
topology.  Electron  micrographs  on  cermet  films  prepared  by  the  same  method  support 
this  conclusion  [7,8]. 
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Room  temperature  relative  reflection  data  were  obtained  using  a  Nicolet  740  Fourier 
transform  infrared  spectrometer  with  a  nitrogen  purge.  The  spectral  range  of  400  to 
1 5000  cm"  was  covered  using  three  combinations  of  sources,  beamsplitters,  and  detectors. 
The  resolution  was  4  cm"  * .  The  samples  were  mounted  on  a  Spectra  Tech  Model  500 
variable  angle  specular  reflectometer.  The  angle  of  incidence  was  5°  from  the  normal. 
A  100%  Pt  film  from  the  set  of  samples  was  used  as  the  reference. 

THEORY 


The  predictions  of  three  effective  medium  theories,  the  Maxwell-Garnett  (MG)  [9), 
Bruggeman  (BR)  [10],  and  the  Sheng  probabilistic  growth  model  (PG)  [3],  were  compared 
with  the  data.  Each  model  is  based  on  a  specific  microstructural  topology.  In  using 
these  models,  we  assume  spherical  grains  and  homogeneity  of  the  cermet  on  the  scale  of 
a  wavelength. 

The  topology  of  the  MG  model  [9]  is  a  dielectric-coated  metal  sphere  (one  may  also 
choose  a  metal-coated  dielectric  sphere).  The  thickness  of  the  coating  is  determined  by 
the  metallic  volume  fraction  of  the  medium,  f.  The  effective  dielectric  function  is 


t\tG  =  f-o' 


2t0(l  -  f)  +  ex(2f  +  1) 
(■o(f  +  2)  -f  £,(1  —  /) 


Sd  r  P  dlel'c,ric  f^ct'ons  of  the  host  and  metal,  respectively, 

is  nn  y  a5ymmetnc  '««*  treatment  of  the  constituent  materials.  There 

no  percolation  transition  since  the  metal  spheres  have  an  insulating  coating  for  any  f. 

dielectric  sphered The  !!",  Cer!"*i  “  3  °f  “ncoated  metal  and 

dielectric  spheres.  The  effective  dielectric  function  is  obtained  by  solving 
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This  model  treats  the  constituents  in  a  symmetric  fashion.  The  percolation  threshold  is 
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where  pj  is  the  fraction  of  oxide-coated  metal  spheres  in  the  composite,  pj  is  related  \> 

to  f  by  means  of  a  free  volume  argument  by 

(l  -  /1/3)3 

Px  “  ,1  -/‘/Y  +  d  '(1  ~/)l/3)3‘  (n 
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<1  and  «2  are  the  dielectric  functions  for  the  oxide-coated  metal  spheres  and  the  metal- 
coated  oxide  spheres  [11],  respectively.  This  model  predicts  percolation  at  fc  =  1/2. 

The  model  treats  the  reflection  of  normally  incident  electromagnetic  radiation  by  a 
thin  film  on  a  substrate  according  to  standard  methods  [12,13].  Multiple  reflections 
within  the  substrate  crystal  were  treated  as  incoherent.  Published  optical  data  on  bulk 
Pt  [14]  and  evaporated  AI2O3  films  [15]  were  used  to  calculate  the  optical  properties  of 
the  cermet  film.  The  substrate  crystal  was  assumed  to  be  aligned  with  the  axis  of 
symmetry  normal  to  the  surface.  The  oscillator  parameters  from  Barker’s  fit  to 
reflectivity  data  [16]  were  used  to  model  the  complex  dielectric  function  of  the  substrate. 

RESULTS  AND  DISCUSSION 


Figure  1  compares  the  measured  and  calculated  relative  front  reflection  for  three 
selected  films.  Figure  la  shows  that  both  the  MG  and  PG  models  provide  good  fits  to 
the  measured  relative  reflectivity  of  a  f=.23  film,  which  has  the  lowest  metallic  volume 
fraction  of  any  film  in  our  set.  The  MG  and  PG  models  make  similar  predictions  at  this 
relatively  low  f  because  most  of  the  coated  spheres  in  the  PG  model  are  dielectric-coated 
metal  spheres.  The  oscillatory  nature  of  the  reflectivity  arises  from  an  interference 
effect  due  to  multiple  reflections  within  the  cermet  film.  These  oscillations  are  less 
apparent  in  films  with  greater  metallic  composition.  To  obtain  the  fit,  the  thickness  of 
the  film  was  assumed  to  be  2200  A,  whereas  the  measured  value  was  1525  A.  The 
structure  below  1500  cm*1  is  due  to  optic  phonons  in  the  AI2O3  host  and  substrate. 

Figure  lb  shows  the  results  for  a  film  with  f=.58.  This  volume  fraction  is  near  the 
percolation  threshold  for  this  system  [6].  The  measured  film  thickness  of  1475  A  was 
used  in  the  calculations.  None  of  the  models  describes  the  data  very  well.  This  result 
is  in  agreement  with  recent  work  showing  that  effective  medium  theories  do  not  apply 
near  percolation  [4],  Details  of  cluster  morphology  not  included  in  the  effective  medium 
treatment,  such  as  nonspherical  shape,  are  most  important  near  fc.  The  rather  flat 
frequency  dependence  of  the  relative  reflectivity  is  also  a  characteristic  feature  near  fc 
[4,17],  The  MG  model  best  fits  the  magnitude  of  relative  reflectivity,  but  the  BR  and  PG 
models  provide  a  better  description  of  the  frequency  dependence. 


Figure  1c  shows  that  the  BR  and  PG  models  provide  the  best  fits  to  a  1300  A  thick 
film  with  f=.83.  The  agreement  with  the  MG  model  is  poor.  A  reverse  MG  model  using 
metal-coated  dielectric  spheres  would  provide  a  better  description.  In  fact,  this  is  the 
dominant  topology  of  the  PG  model  for  f  near  unity. 

Berthier  and  Lafait  [18]  have  studied  the  infrared  properties  of  Pt/Al203  cermet 
films  prepared  by  rf  cosputtering.  They  were  able  to  fit  the  measured  optical  constants 
using  the  BR  model  over  a  broad  range  of  composition.  However,  they  treated  the 
volume  fraction  as  an  adjustable  parameter  and  introduced  an  effective  depolarization 
factor  as  well.  The  BR  model,  with  fewer  adjustable  parameters,  does  not  apply  to  our 
films,  which  were  prepared  by  a  different  method. 

In  conclusion,  we  find  that  the  PG  model  provides  the  best  description,  among  three 
simple  effective  medium  theories,  of  the  relative  reflectivity  of  thin  Pt/Al203  granular 
cermet  films  over  the  complete  range  of  composition,  although  the  agreement  is  not 
quantitative,  particularly  near  fc.  We  intend  to  extend  this  work  by  measuring  the 
absolute  reflectivity  using  the  W-V  method  and  to  extract  the  frequency  dependent  ,.j 

complex  index  of  refraction  of  films  for  which  the  transmission  is  also  measurable. 
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